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A VISUAL STUDY OF RADIAL INWARD CHOKED FLOW 
OF LIQUID NITROGEN 

by R. C. Hendricks, R. J. Simoneau, Y. Y. Hsu 
Lewis Research Center 


ABSTRACT 

Data and high speed movies were acquired on pressurized 
subcooled liquid nitrogen flowing radially inward through a 

0.0076 cm . (3-iiiil) gap. The stagnation pressure ranged from 
0.7 to 4MN/m (0.2 <_ P/.P £ 1.2). The ..results of this quali- ■ 
tative study indicate: 

1. Steady radial inward choked flow appears ' equivalent 
to steady choked flow through axisymmetric nozzles. 

2. ‘- Transient choked flows through the radial gap are 
not uniform and the discharge pattern appears as nonuniform 
impinging jets. 

3. The critical mass flow rate data for the transient 
case. appear different from those for the steady case. On the 
mass flow rate vs pressure map, the slope and separation of 
the isotherms appear to be less for transient than for steady-. 
radial choked flow. 


SYMBOL LIST 

- 2 •■■■' . • d' ' . ■ ■ 

area , cm ... 

hydraulic diameter, cm 
mass flux, g/ cm^ -sec 
reduced mass flux 

2 

mass flux normalizing' constant g/cm -sec. 
For nitrogen G*=6010. 

passage length, cm 

mass , g 

pressure I-JN/m (psia) 


A ' 
D ' 
G 



L 

M 

P 
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P p =P/P c reduced, pressure 

T temperature, K 

Tp=T/T c reduced temperature 

3 

V volume, cm . 

y vertical distance, cm 

Z=P/pRT compressibility factor 
A . difference 

p density, g/cm 

r time, sec 

Subscripts 
bulk . bulk 

c critical 

o stagnation 

V-' 

sat saturation 

R reduced 


INTRODUCTION 

Choked flows are common in gas lubricated bearings , seals 
and in the venting of high energy systems when they "go 
critical." - However, little ..is. known ..about using a vaporizing 
liquid as a seal or in bearings, or leaks through cracks in 
either the transient or steady state cases. 

The motivating problem of interest is to be able to pre- 
dict the leak rate and possible onset of two -phase flow in- 
stability due to maldistribution of pressurized cryogens dis- 
charging into space through a long narrow rotating passa.ge 
with.heat being added along the walls and/or generated by the 
fluid. 

Some of the unique features of the problem are: the small 
gap heights required CO* 0008-0 .0013 . cm) , L/D values of 
(100-200) , and a flow aspect ratio (passage front/passage 



height) of order 10 . These parameters could result in mal- 
distribution of vanor and liquid flow, surface capillary 
effects, thermodynamic nonequilibrium, boiling and vapor 
entrainment, etc. While the efforts of several researchers 
have been directed .at two-phase choked, flows (refs. i } 2, 3, 
and 4). Little or no research has. been done on these ques- 
tions within the above range of parameters. 

A quantitative study was undertaken to delineate some of 
the problems associated .with: radial inward flow through a 
narrow stationary gap which, chokes, and discharges through a 
central pipe. The data and accompanying high speed movie 
were, taken to explore the problems. 


APPARATUS 

The basic apparatus used in this study is illustrated 
in figure 1. Pressurized liquid nitrogen flowed through the 
radial gap and out of the cryostat. In this study, the 
radial flow passage was formed by a glass cover separated 
from the stainless steel orifice by three radial wires 
0.0076 cm diameter £3 mil). Details of the radial gap-nozzle 
are given on figure 1. . • 

A 10 -common-wire Chrome! -Alumel -.thermocouple r.ake. was . . 
constructed to monitor ..the -liquid level in the reservoir.. 

The rake frame was made of, bake life and ..the thermocouple 
junctions were made at 0.63 cm (1/ 4 in.) intervals from the 
face of the bakelite insert. 

The reservoir temperature was monitored by a platinum 
thermometer, located 0.. 95 cm £3/8 in.) above the radial gap. 

A second platinum thermometer near- the top of the thermo- . 
couple rake was used as a check under no flow conditions. 

The vacuum jacket was usually not -better than 50 Op, how- 
ever, some of the beat leak was offset by the evaporated 
nitrogen which leaked past the window seals. 

STEADY -STATE RESULTS 

.In order to obtain steady-state critical flow data, the 
entire apparatus .of .figure 1 was mounted in the blowdown- 
facility’ (f.ig,. 2) with a reservoir capacity of 120 liters 
which is quite sufficient, for all steady flows taken, in this 
study. . The entire system could be precoo.led. from ambient 
temperature down to about 80K, which enabled us to achieve 
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various reduced temperatures starting at. T f , '= 0.7. Since the 
blowdown facility: was. sized, for- an order of magnitude larger 
flow area, the extremely: low: mass, flows could not be metered 
to a certainty greater than 5 percent; however:, for. gaseous 
flows the two flowmeters, agreed within 3.. percent. Peculiar 
temperature inversions were often noted in the flows from the 
reservoir of figure 2 into the apparatus of figure' 1. These 
may have been due to tank •stratification, coupled with gas 
residual in the. ..apparatus „ For.: all the data reported herein, 
.the flowmeter between the reservoir and the apparatus was used. 

The critical, mass flow rates . for several fluids can be 
correlated using .. the principle of corresponding states , ref er- 
ence. 4. The flow rate normalizing parameter 

G * - c<-cW 1/2 v 

reference 4, is used to normalize .the. flow data presented 
herein. The. gaseous .data, figure 3 ; and table I, appear to be 
in reasonable agreement with ..previously determined data for 
gaseous flow through nozzles. This, of course, gives us a 
certain degree of confidence .-in our -results . The trends in 
the data for the T. R .= 0..8 isotherm are very similar to those 
.of the axisymmetric nozzle data (refs. 3. and 4). Solid lines 
on figure 3 are used. to represent the data of reference 4; 

Data trends for other isotherms are in equally good' agreement, 
as can be noted in figure. 3. One questionable' area appears to 
be the TV = 1.035 isotherm near = 1.2.' These data are 
presented in table I. 

The general agreement. between the data of table I and 
that. of references 3 and 4 indicate that there, is little 
difference in, the' critical flow phenomenon between the axi- ■ 
symmetric nozzle and radial inward flow through a 0.0076 cm 
(3 mil) gap having an outer, to .inner area ratio of. «* 13:4,' 
see figure 1. This is quite significant as the geometries are 
very different and tends to imply the generality or the useful- 
ness, in the application of the technique of reference 4 to. odd 
geometries. 


TRANSIENT RESULTS 

Because of . the ..uncertainties, in determining the parameters 
which make ...lip. critical, f low rate , i . e. density, area, and 
.volumetric flow, rate, exact values of G were not determined.' 
However some interpretations of the data were made and are 
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presented in the appendix. In essence,, two . techniques- are 
used to predict transient data with limited success. - The 
transient results depart in', both slope , arid level when Com-, 
pared to the . steady, state results. Further investigation 
of the disparities are required. 

PHOTOGRAPHIC SEQUENCES 
Transient Conditions 


Motion pictures were taken at 100, 400, 5000,' and 16 200 
frames per second to provide greater detail to flow patterns 
and the gas-liquid interface within the reservoir. For ail 
tests, the back pressure is near ambient and the fluid is 
nitrogen. 


P critical = 3 - 417 MN/m C33.7.2 *til0 ; 

T • . • , = 126.3 K; 

critical , ’ 

»ctittcai --'0.3105 -gm/ce. , 


During the flow transient, the ullage space is continually 

charged with high .pressure, nitrogen ,gas. at ambient . tempera- 
ture. The charge fate is not: sufficient ..to maintain cons tant 
pressure but does provide adequate isobaric sUbcobling* i.e. 

^s at 4 bulk^ > V . 

The combination of decreasing- pressure, with increasing temper 
atiir.e provides. ..a spectrum, of: f.lo.w regimes from subcooled , to 
saturated i to gas-dropiet entrainment. 

in general . the.-.- f low patterns are irregular and the f low 
area covered, by the single . phase, fluid is. much greater than 
the area where two -phase phenomena occur. 


Description of Setup • 

Sequences to illustrate the., .behavior of - the gas- liquid 
interface, during transient flow were taken at 100 frames per 
second. The camera viewed the,, interf ace through the side, 
windows of the cryostat, and the apparatus is as ' illustrated 
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In figure 1. The thermocouple rake , , platinum thermometers , 
glass cover plate as well as. reservoir thermals are clearly 
illustrated.- These sequences could not be taken simultane- 
ously with those viewed from the top. 

To obtain. more definitive pictures of radial inward . 
choked flow as viewed from the 'top and raise the liquid level, 
the single 3/8 in. -thick glass cover plate of figure 1 was 
replaced by a stack of three 3/8 in. -thick glass plates, and 
the. camera .speed increased to 400 frames per second. . The 
pictures were taken while maintaining . the .upper portion of 
the tank.,, namely, the ..camera .viewing ...windoxv, above the ..conden- 
sation temperature, and .filling the dewar only to the level 
of , the.; stacked...glas.s . plates . . Previous, experience indicated 
that., ullage gas would condense on the window and other., parts 
of the . dewar : form liquid drops and/or rivulets which obscured 
direct observation. 

In the quest. for. more definition, top-viewed sequences 
were taken xvrth... top and side lighting at an effective, camera 
speed of 16 200 . frames ..per. ..second with an. equivalent shutter 
speed of 152 000 . frames per second. This x^as achieved, by 
running the. camera at .81.00 f rames..._per ,s.e.cond and recording 
an. 8 mm -s.pl it frame.. ...image ...as ...viewed. ./through, a-, narrow- slit. 

At these speeds . the entrainment of. voids can: be studied but 
camera speed of another, order of magnitude is required to 
study the choking interface. 


Description, of Film Sequences 


. Film Sequence 1 

To illustrate the flow patterns in. s.ubcooled-subcritical- 
nitrogen the- stagnation., pressure and; temperature in the. cryo- 
stat . is . set at. 0,.. 7.9. MN/m 2 . (114. 2 .psia) and 87.1 K respectively 
which represents . a .subcooling of 13 K. At the opening of the 
valve:, the .fluid surf ace vibrates, and. a. white, central two- 
phase; region begins-., to. form. Mote that the f loxv. discharges, 
as an 'array .of .nonunif.orm. lets .. ....,The ..f low patterns are not 
stable, and ..the . jets, extend, inward far. beyond: the nozzle edge, 
the. minimum area., ...where choking. ..should, occur... Such behavior 
.indicates a substantial ...degree of...noneaui librium. The average 
pressure during., discharge . .is .-. 0....-7 5 . MN/mr. (10.9 ps.ia) with a 
depressurization rate of -0.0035 MM/mr-sec ( - 0.5 psia/sec) ... 
While the .*g.'as,r,liq.ui.d..,interf.ace.. remains . quite calm, thermals 
as. well as. .flow., patterns of., the., fluid approaching the. nozzle 
can be seen. At the reduced- liquid levels, the temperature of 
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the f luid entering the, nozzle approaches, .saturation and the 
. .jet sources .begin . to . extend ,.o.u twar.d. The.. ; effects .. of. small' 

.. s.ur.face;,defe.c.ts .', Gp its. .or„sxte, si)_can..be\. noted., in, .addition .to 

. the nonuniformity of.:. the. .two -phase front,... Further .reduction 

in tlxe level permits, .the,, visualization ..of strong surface, 
tension., ef fects.,... the . entrainment.. ef,. gas..-. and ...a , .rapid develop- 
ment of a fan pattern which finally, fills the radial gap- 
nozzle. 

Figure 4(a) represents this film sequence. 


Film Sequence 2 

In direct contrast to the subcrltical pressure patterns 
.of Film Sequence 1, .this... sequence represents., fluid ..nitrogen 
in. its _ supercritical., pressure state . . The stagnation pressure 

and temperature are 4 .2 Ml/mr .. 16.0.9,.. 2..,ps.ia3 . and. .121. 8 K 

respectively. Here. .the -.transposed , critical temperature* is 
130 . 7 K, so the fluid ,1s..., in. .a sense sub cooled nearly. 9 K... At 

the opening of the. valve.., . the ...central two-phase region (white) 

develops rapidly. The region, is ..not . steady nor does it have 

.periodic, oscillations and distinct ..j.etting is difficult to 
determine . The average pressure. during discharge is 

4.0 MN./m. 2 f 5 7.9 psia), with... a ..depressurization., rate of 

- 0 . 4 7 MN /m 2 -sec . (- 6 8 ps.i a/sec) .during . the., ini t i al . phase of 

the. transient .and. „rd,.,0.9.6.,MM/m 2 . -..sec, ^C--:l 3. 9 psia/ sec.)., during 

the ma.j.or ..portion ,.o.f....f he-,.discharge..':cy.c.le . ,:,. A..., view, of the gas- 

fluid. interface.. reveals the format i onlo.f ... .a I'm i.l ky.!. ? layer 

. approximately 1 cm, -thick.... B.y . carefully..! ollowing the proce- 

- dure, described, in . Descr.ip.tion of Setup., the th ickness .. of this 
.layer was minimized and pictures of the flow patterns were 
obtained. ' . • 

Figure’ 4 Cb.) represents:, this.. .f ilm-sequence... • As the. run 
progresses , . the. .reservoir pressure, decreases and the liquid 
temperature Increases due . to heat,, transf er., through the.-, ullage 

gas,, from the. make... up, gas.. At. the..,beginning. of, -the fan. pattern 

the. pressure., is 3 . 86 MN/.m“- ., (5.4.5,. 2. psia) ..and. the' temperature 
measured, a.t., the. .platinum, thermometer, level which is; 0...95 cm. 

(3/8 in. ) above the, radial flow.. plane, . is 129. 1. K. This temper- 

... ature ,is, .characteristic .... of .. the. mixed , gas -condensate interface ■ 

... . . .... *The, transposed, critical temperature, is, a . point corres - 
, ponding. :.to, a... sharp, p.eak.’.in,.s.pe.cif.lc.: heat.,,. C_..,. . It. is a function 
. of, pressure... At .this. ..temperature the. density, changes ..rapidly 
from. high. to. low : density. . . In. .a. sense, it can. be thought of as 
a "supercritical saturation temperature equivalent." 



and the fan. patterns. which- follow are. presumed to be "mist 
flow.:" Note the multiplicity of variations in the .choking 
front as temperature increases and. pressure, decreases. ; The 
appearance . of inner and outer ring patterns . suggest that two 
fronts exist in the nozzle. At least the picture suggests 
multiple regions with large gradients. 


Sequence 5 

2 

During the course of a run at the 4 MN/m pressure level, 
a nonhomogenity which looked like a bubble, was entrained into 
the nozzle. The reduced pressure and temperature are approxi- 
mately 1.2 and 0.96 respectively. Bubbles should not form, 
under these conditions. The picture sequence is expanded in 
order to view the nonhomogenity. . 

. .... Sequence 4 

. To illustrate flow. patterns near the thermodynamic 
critical point, the initial pressure is 3. 4 8 l!N/m, (506.2 
psia) and the temperature is. 126.9 K. These conditions are 
just above the thermodynamic. critical point. The temperature 
is quite close to the transposed . critical temperature of 127 K, 
consequently the "subcooling" is smaller than the previous 
sequence. Note the aperiodic oscillations of the two-phase 
region which appears white. Further note that on an area basis, 
most of the flow field within. the nozzle appears nonoscillatory 
and homogeneous. The average pressure during discharge is .7 . 
3. 3 .MN/irw (480 psia) with an initial depressurization rate of 
-0 . 3:8. MN/irw-sec (-54.5 psia/sec) and average of -6. 092 MN/nw-sec 
(-13. 4 psia/sec) . As the reservoir empties the fan' pattern is 
observed. The pressure has decreased to 3.26 \MN./irw (472.2, 
psia) .. which is less than critical pressure while' the temphr- ' 
ature has increased.. to. 127.8 IC.at the platinum thermometer, 

0.95 cm above the nozzle plane. The mist flow regimes become 
apparent, even though the conditions are quite,, close to the 
critical point. 


Sequence 5 ' ' 

In the following three sequences the flow patterns are 
illustrated for three selected subcooled-subcritical pressures. 
At the opening of the valve the pressure is 2.82. 1 HI (409.2 • 

psia) and the temperature is 113.6 K. This represents a sub- 
cooling of . 9 . 4 K. This level of subcooling persists over most 
of the film sequence, however at the development of the fan. 
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pattern tlie temperature . is . within.fi ..K of the saturation 
temperature. The. average -pressure is 2 . 7-..MN./mr.-' C 3 S 4 psia) 
with an initial depressurization. rate of -0.22 MN/m. -sec 
(-32.5 psia/sec). and an average rate of -0. 045 MN/m^-sec 
(-6.5 psia/sec) . The oscillations of the two-phase region 
become more apparent and. the choking front begins to extend 
outward, however the majority of. the nozzle is still in 
homogeneous nonoscillating flow. At the onset of the fan 
pattern evidence of an interface or. surface tension effects 
on the entrained flow can be detected. 

As time progresses, the entrained mist flow transitions 
to high density vapor flow. 


Sequence 6 . 

2 

As the valve opens the pressure is 2.16 MN/m (313.2 
psia) a:nd the temperature is 96. 8 It.. This., .represents, a sub- 
cooling of 10.2 K. As : the flow rate has diminished consider- 
ably, the ^pressure-temperature- time profiles are nearly 
linear. The average pressure is 2- MN/m° (294 psia) and the 
depressurization rate is -0 . 037 MN/m 2 - sec (-5.3 psia/sec). 
Oscillation in the two-phase region is commonplace and the 
nonuniform two -phase front ..has . extended in an erratic manner 
into the. nozzle. Most of the nozzle on an area basis is still 
in uniform, flow. As the radial fans appear, the effects of 
surface tension. are now easily observed. The pressure has 
decreased to 1.95 MN/nw (283.2 psia) and the temperature at 
the platinum thermometer above the nozzle plane is 117.5 K. 
Again the nozzle flow transitions through several flow re- 
gimes from subcooled saturation, super saturation, mist flow 
and finally cold vapor.. Wi thin, these latter frames- one can 
visualize . what the nozzle flow would look like., if it . were 
flowing-- steady state--in each of these regimes furthermore 
the picture sequence illustrates the transition boundaries 
between these states. Peripheral waviness noted in the. pic- 
tures is due to light ref lections .from the interface and does 
not represent nozzle characteristics. - , 


Sequence 7 

At 1.5 MN/m 2 (217.2 psia) and 9 4.3 It the . reservoir is 
subcooled 5.9 K. As the valve. opens, the development of the ^ 
choking region, is observed.. The average pressure is l.^MN/im 
(205 psia) .with a depressurization rate of -0.0147 MN/m -sec 
(-2.1*2 psia/sec). The two -phase region is oscillatory in 
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nature and the front extends erratically into the nozzle. 

While most of the nozzle appears homogeneous., the erratic two- 
phase radial spikes seem to make . the .uniform circular , two- 
phase choking front .model a mathematical concept whose useful- 
ness is fortuitous. The time temperature -nressure profiles 
are quite linear and diminish. to 1.39 MN/m^ (201.2 psia) and 
109 K at the . platinum. above . the. nozzle .plane as . the fan pat- 
terns begin. The effects of surface tension are clear and 
appear to exert considerable forces on. the flow. field. The 
nozzle. is clearly not in equilibrium . and. combinations of drop- 
lets, cold vapor and liquid are seen to. coexist within various 
sections. 


Sequence 8 

To illustrate the flow characteristics of a saturated 
fluid, the 1.5 .MIT./ m ^ pressure was chosen because of its rela- 
tion to a seals project at Lewis. Research Center. To view the 
phenomena more clearly' the frame .rate was /changed, to about. 

5000 frames per second and the nozzle lighted from the top and 
the side. Mote the extreme irregularity of the interface 
which is delineated by the. two-phase streaks toward, the nozzle 
center. It would seem ..that the . two -phase choking interface 
per se does not exist but rather a highly nonequilibrium non- 
steady transition between . liquid., and vapor . As bubbles are 
ingested,, large vapor rays separate the liquid and plunge into 
the orifice exit. • 


Sequence 9 

In order to study these bubble streaks... -better , the frame, 
rate is increased to 8100 8 mm-split. frame images per second, 
equivalent to 16 .200 £rames/sec. .One can. now. see .that: bubbles 
are elongated. as they approach the entrance of . the nozzle and 
are sucked through. The leader seems to be spherical and non- 
steady, but progresses in a. linear manner to the "interface" 
or choked front. In this stop frame sequence one can see the 
advance of a bubble about 30° off. the. bottom. of the picture 
(counter- clockwise) . and. shortly after this bubble is halfway 
across, the nozzle a . second bubble enters at 10° (clockx\rise) 
from the bottom. The vapor streamers persist for quite some 
time before the gap is closed. 

Figure. 4(c) represents this film sequence. 



Sequence 10 

, To illustrate the...inf luence of .surface, pits and. scratches, 
we see here the activation, of . a .si.te . about: ;30° counter- cloclavis 
from the bottom..: and .midday across the nozzle. The camera speed 
is 16 200 frames per second. ' ■ , ■ 


SUMMARY -MOTION PICTURES 

■In- summary., motion pictures have been';.taken .which illus- 
trate 1 several aspects of radial inward choked flow of fluid 
nitrogen. . 


Sequence 11 

This split frame sequence compares radial inward, choked ■. 
flow at supercritical and subcri.tical pressures.; The: reduced 
pressures are approximately P^ = 1.17 and Pjv = 0.22 respec- 
tively. It does appear that . the flow field^is more homogeneous 
for P R >1; because of' its erratic nature the "choking!', inter- 
face, appears larger for .Pp ; < 1. Both flow fields are. oscil- 
latory; for P jV < 1 the., amplitude . and oscillations are more 
clearly -defined. . • • • • 

' • . Sequence T2 

The erratic nature of two -phase- interfaces leads one to 
believe .that choking. . is . a . result of nonequilibrium nonsteady 
transition from liquid to vapor. 


Sequence 15 

Motion pictures at 16 200 frames/sec reveal that bubbles 
are elongated and . sucked.. through., the gap in a radial manner.. 
The leader appears spherical: and. nonsteady. . Higher frame 
rates are necessary to study the choked interface. 


Sequence 14 

A wide variety of flow regimes (subcooled, saturated, 
separated, mist, and transitions between these regimes) were 
illustrated as the liquid level dropped. 
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CONCLUSIONS 

The quanti tative .. study undertaken . to determine some of 
the . problems;, associated x^ith discharging pressurized cryogens 
through a radial' gap has revealed that: 

1. • Steady radial inx^ard critical mass flow rates appear • 

equivalent to critical mass flow rates through axisymmetric 
nozzles. ' 

2. . The radial flow of two-phase and supercritical mix- 

tures through, a narrow gap is no.t uniform. For subcooled 
flows the- surface area covered by single phase. fluid. is . much . 
greater than that covered by the two-phase discharge. As • sub- 
cool ing diminishes the tx\ro-phase area .approaches the single 
phase area. Over a broad range. in quality (including the near 
critical region), the basic flow patterns are made of radial 
jet streams which may occupy a segment of . the radial gap- 
nozzle or the whole gap area. Flow maldistribution is quite 
common and oscillations in circumferential and radial direc- 
tion are obvious;. ‘ , • • . . • . • 

i'-. > 4 • 

3. . The critical mass flow rate data- for the transient 
case appear different from those' for the steady case.- On 
the mass -flow rate as a function, of pressure map. the. slope 
and separation of isotherms for transient radial choked floxv , 
appear, to be less than for steady radial choked flow.’ Further 
work and quantitative data are required to resolve the dis- 
parities between steady and transient critical radial inward 
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APPENDIX-TRANSIENT RESULTS 

The reservoir, temperature (lower thermometer of fig.,1) 
pressure and thermocouple rake were monitored .on . strip charts. 
To record data, the reservoir was filled and .pressurized; 
the recorders started, and an . electronic, timing marl; placed 
on each chart; the exit valve was then opened. At low re- . 
duced pressures, a loud "crack" signaled. the onset of flow 
out of-, the vent stack. The liquid . level could be monitored 
visually (even for P > P. ) and the. run was terminated when 
the liquicl-gas inter.f ace reached the lower part of the glass, 
see figure 1. . 

The thermocouple rake, .was calibrated using 1-1-1 tri- 
chloroethane and observing the merging. ..of., the refracted, and 
reflected images. of the .wire while metering . the expelled, 
liquid. With this volumetric calibration, and . the response 
of the thermocouple rake, the. fluid volume. (V) exhausted as. . 
a function of time . (x) was determined. Evaporation of the 
liquid off the thermocouple wire as the . liquid vapor inter- 
face passed usually gave a sharply . defined constant time 

temperature trace from which one could . determine the dewar 

level variation with time.,, see figure 5. Either the onset 
or termination of evaporation liquid and gas. indicators 
respectively, could be used . to .determine the ..volume (V) vs 
time curves; generally, the onset .of ... evaporation, •(liquid 
indicator) was. used. A typical series of dewar. volume vs 
time, curves are given. .in f igure 6 ... Although the ■ data do not 
always indicate a linear ..relation between Ah and x , the uncer- 
tainties did not .merit the use of a higher degree curve. It 
appears that for a given 

C?o> T o).>. = constant , . .. ' . ' ;• ' • ; y 

gives a nearly accurate description of the volumetric change 1 
in the dewar. •’ 

From. the experimental data 


(? 



J 



values for the transient critical mass . flow rates can be 
calculated and compared. with those for. the steady flow case. 
If the critical mass flow rate is expressed as 


G 


R 


C-dM/d-O/AC-O 

G* 


G 

G^ 


max 


CD 
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where. 14 = pV, the mass of, fluid remaining in the reservoir 
and A(t) is the area of the choking front one can - see immedi- 
ately some of .the problems in interpreting the data. ; Even 
though the rate of change of volume was reasonably constant , 
the average density of the fluid is difficult to determine 
because the pressure and temperature change with time, and 
temperature also may change in space. A value for the average 
density was calculated assuming that the temperature-time 
profile also represented the temperature-space profiled 



The density; corresponding to the stagnation pressure, and . 
assumed temperature was found using a generalized thermo- .; 
physical properties program GASP (ref. 5). Typical histories 
for reduced stagnation pressure and temperature are given, in 
figure 7. Note that this procedure assumes the bulk fluid to 
be in equilibrium, i.e. density uniquely determined by the . 
instantaneous values of pressure' and temperature in the bulk; 
an assumption not necessarily valid in the transient. case. 
Finally, the actual location and area of the choking front 
cannot be instantaneously determined. Consequently the abso- 
lute vdlues for transient critical mass flow rate cannot be 
defended;, They do provide the basis for the comparison with 
tile steady flow case, where again the actual choking front 
area had to be assumed. The value used for A was based on 
the inner diameter and passage height. 1 . 

A. typical variation of G/G* = Gp with time is illustrated 
on figure 7. Values of G/G* . =. Gr along . selected isotherms 
are shown ori figure 8, as a function of reduced . pressure . 

One must note that these profiles, (isotherms) appear to be 
nearly independent of changes with Tr, while the isotherms; . . 
of figure 3 show definite changes with. Tr.- Furthermore, . the 
slope of these isotherms (AGr/APr) is considerably less than 
for the steady-state case; figure 3. 

An, even simpler approximation to the transient, reduced 
critical flow rate is shown, in. figure 9. Here the fluid, 
.'density at any time is simpiy taken to be that calculated 
from the, stagnation pressure and temperature, f 0 (t) and P 0 ’(t) , 
rather than attempting to correct for spatial variations 
through the fluid volume. 
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A typical variation o£ equation. (3) with' time, is given 
in figure 7. While .ithe slopes , of the isotherms on figures .: 
3... and, -9 are nearly the. same , the level .change... between iso- . 
therms is. much, greater;. In. figure 10.,- a comparison is made, 
between the. steady results , f igure 3 , and the .'transient re- 
sults, figure 9 ,, for reduced temperatures of O.'o. and .0 . 9 . 
Deviations in both slope and levels are to. be noted. 


The discrepancies. .may . arise ..from.. .the diagnostic limi- 
tations of the, present , .basically qualitative experiment . 

They, may on. .the other hand stem' from phenomena associated 
with the transient conditionr-such as nonequilibrium density 
in the fluid, or rapid changes in location and area of the 
choking front. In either. case, further .experimentation with • 
more elaborate apparatus is needed to resolve the question. 

, The transient data obtained in this investigation are 
presented in table IT.- . ■ . 
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" 'TABLE L. -RADIAL ■ I MW APvD. C RI T IC A L FLOW 
' STEADY STATE DATA 


Run 

P R 


T 

X R 

g r 

407 

0.60 3 


2.24 

. ' 0. 0909 

40 8 

.603 

- ■ ; ; 

2.24 

..0913 

409 

1.202 


. 705 

a l . 18 

410 

1.022 


.706 

• : ■ 1.0 56 

411 

. 816 


.70 8 

- • . . . 907 

412 

. 606 


.712 

'. 766 

413 

. 426 


.717 

.615 

414 

. 255 


.718 

.377 

415 

.2.51 


.719 

. 391 . 

416 

1.218 

r 

.793 

1.107 

417 

1.032 


.793 

.942 

418 

. 818 


.798 

.814 

419 

. 611 


. 799 

. 666 , 

420 

.411 


. 805 

! ; 423 

421 

. .248 


. .-804 

. 0342 

423 ‘ 

1.206 


.918 

. 818 

424 

1.019 ■ 


. 898 

. 729 

425 • 

.815 


.900 

. 603 

426 

. 59 


.902 

. 391 

427 

.572 


.900 

.140 

428 

'.556 


. 891 

a .062 

429 

.485 


. 891 

. 058 

430 

-.599 


2.0 5 

.0907 

431 

.599 


2.05 

. 0904 

432 

1. 218 


. 867 

.950 

433 

.985 


. 869 

.800, 

434 

-.794 


• .870 

- : .659 

435 

. 592 


. 874 

.475 

436 

. 592 


. 872 

. 475 

437 

.445 


. 877 

. .0945 

438 

1.190 


.988 

. 568 

439 

1. 121 


. 96 7 

. 589 

440 

1.117 ' 


.966 

. 5895 

441 

1.034 


.965 

. 523 

442 

.995 


.965 

. .489 

443 

. 958 ' 


.966 

. 459 


Questionable. • 
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TABLE I. -CONCLUDED. -RADIAL INWARD CRITICAL FLOW 
' ■ STEADY STATE DATA 


Run ' 

P R • 

T 

1 Pv 

g r 

444 

0.878 

0.970 

0.192 

445 

. 1.266 

.9 75 

a . 76 

446 

• 1.267 

1.044 

. 314 

447 

1.267 

1.042 

. 314 

448 

1.205 

1.032 

.291 

449 

1.128 : 

1.023 

. 262 

450 

* 1.052 

1.011 

.236 

451 

1.025 

1.007 

.226 

452 

.927 

. 990 

.190 

453 

. 844 

.975 

. 163 

454 

.646 

.934 

’ .101 


Questionable . 



TABLE II. -TRANSIENT RADIAL INWARD CRITICAL FLOW 
DATA INTERPOLATED TO SELECTED ISOTHERMS 
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TABLE II. -CONCLUDED. -TRANSIENT. RADIAL INWARD CRITICAL FLOW 
DATA INTERPOLATED TO SELECTED ISOTHERMS. 


20 







tO tO 03. VO 

to CN1 

VO rH 

LO CD 





cm 

VO [X* 

VO 

l--* 

IX* 

00 


CO 

CO 





pi 

• • 

• 

• 

• 

• 

• 

• 

• • 






O 













to r-l 

CD 

tx* 

(X] 

to 

LO 


CM CD 





rH 

CO cd 

CO 

CD 

CD 

CD 

CD 

rH O O 





pi 

• • 


• 

• 

• 

. 

• 

• • 




■ 

CD 

o 






rH 

rH rH 




cn 













• 

H 

(XI 

to 

CO 

1 

6 

2 

4 

2 

6 




o 

■ <3 

• • 

• 

• 

• 

• 

• 

• 

• • 





„ 

co cd 

LO Xj - 

O 

rH 

LO 

(XI 

CO CD 




n 

> 

CD O 

o 

rH 

rH 

(XJ 

rH 

CM 

CM t-H 





<1 

t-H 

rH 

rH 

rH 

rH 

rH 

rH 

rH rH 




ci 













H 















t-H 


rH 

LO 

LO 'd- *v|- 

LO LO 





Pi' . 

CD CD 

CD 

CD 

CD 

CD 

CD CD CD CD 





CL 

• • 

« 

• 

• 

• 

• 


• • • 






t-H i— 1 

t-H 

rH 

rH 

rH 

rH 

rH 

rH rH 






to rH 

(XI 

CO 


to VO CD 

CM 





Pi 

IX* CO 

CO 

CO 

CD CD CD 

CD 

O O 





p -, 

• . • 

• 

» 

• 

• 

• 

• 

• • 






o 







rH rH 






00 CD 


K 

LO 

CO 

(XI 

CM 

ix* to 


rH 

VO 


CM 

VO VO 

f"- 

t"* 


t"- 

CO 

CO 

CO CO 

00 

CD 

CO 


ci 

• • 

• 

• 

• 

• 

• 

• 

• • 

• 

• 

• 

■ .v 

CD 

o 














to 

LO 

VO 

xj- 


to 

VO LO 


LO vo 


rH 

co co 

co 

o 

CO 

CD 

o 

O CD CD 

rH 

o 

O 


pi 

• • 

• 

• 

• 

• 

• 

• 

• . 

• 

. • 

' • 


CD 

o . 


rH 



rH 

rH 


rH 

rH 

rH 

LO 













CO 













• 


IX* CD 

CM 


to 

IX- 

CO 

LO VO CM 


CM VO 

o • 

< 

• . 

• 

• 

. 

• 

• 

• 

• • 

• 

• 

• 


X, 

to r- 

CO 

LO 

LO 

CD H- 

rH 

rH LO 

CM 

CO 

CD 

II 

> 

! OVCD 

(D 

o 

o o 

rH 

rH CM rH 

(M 

CM 

rH 


<3 



rH 

rH 

rH 

rH 

rH 

rH t-H 

rH 

rH 

rH 

Pi 













H 















CO 00 

cd 


CD 




CM 



(M 


Pi 

o o 

o 

rH 

O 

rH 

rH 

rH 

rH rH 

rH 

rH 

rH 


Q_ 

• • • 

• 

• 

• 

• 

• 

• 

• • 

• 

• 

• 



(XI CM 

CM 

CM 

CM 

CXI 

CM 

CM 

CM CM 

CM 

CM 

CM 



VO LD 

to 

tO CM 

(XJ 

CO 

CM 

LO VO 


CM 

to 


ci 

LO VO 



O 0 

CO 

CO 

CD 

CD CD 

O 

O 

O 


p* 

• • 

• 

• 

• 

• 

• 

• 

• • 

• 

• 

■ • 



o. 








rH 

rH 

rH 



E-7605 




NITROGEN GAS 



I 


n HEAT EXCHANGER 
ORIFICE FLOWMETER 

Figure 2. 


CS-67136 



5 


1.0 

C R 

• .5 

0 



-7605 


Ul 



(c) REPRESENTS FILM SEQUENCE 9. 

Figure 4. - Photographs which represent selected film sequences of the accompanying motion picture. 
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REDUCED PARAMETER. 



T, SEC 


Figure 7. - Variation of parameters with time for an 
initial stagnation pressure of 600 psig. cs . t 
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